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00 ■ ABSTRACT 

\ As part of the European Supernova Collaboration we obtained extensive photome- 

try and spectroscopy of the type la SN 2002dj covering epochs from 11 days before 
lO ' to nearly two years after maximum. Detailed optical and near-infrared observations 

\ show that this object belongs to the class of the high-velocity gradient events as indi- 

OO , cated by Si, S and Ca lines. The light curve shape and velocity evolution of SN 2002dj 

appear to be nearly identical to SN 2002bo. The only significant difference is observed 
in the optical to near-IR colours and a reduced spectral emission beyond 6500 A. For 
high- velocity gradient Type la supernovae, we tentatively identify a faster rise to max- 
imum, a more pronounced inflection in the V and R light curves after maximum and 
a brighter, slower declining late-time B light curve as common photometric properties 
of this class of objects. They also seem to be characterized by a different colour and 
colour evolution with respect to "normal" SNe la. The usual light curve shape param- 
eters do not distinguish these events. Stronger, more blueshifted absorption features 
of intermediate-mass elements and lower temperatures are the most prominent spec- 
troscopic features of Type la supernovae displaying high velocity gradients. It appears 
that these events burn more intermediate-mass elements in the outer layers. Possible 
connections to the metallicity of the progenitor star are explored. 
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1 INTRODUCTION 

In the last decade type la Supernovae (SNe la) have 
been extensively used for cosmology yielding evidence 
that we live in an accelerating Universe (Riess et al. 
1998; Perlmutt er et al. 1999) . Ongoing surveys such 
as ESSENCE (|Miknaitis et al.l l2007l : IWood-Vasev et all 



l2007f ) and SNLS (|Astier et all l2006l 'l use the relation be- 
tween t he shap e of t h e light curve and it s peak lumi- 
nosity jP hillips' '1993'; 'Phillips et al.' '1999*; 'Hamuy et alj 
19961: IRies s. Press & Kirshncr 1996; Goldhabcr et al. 200l 
Guv et al.lf2005l : |Prieto. Rest fc Suntzefl 1120061 ') to constrain 
the equation-of-state parameter for dark energy. An impor- 
tant caveat is that those relations assume that SNe la are a 
one parameter family. 
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However, in recent years there has been growing evi- 
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dence for the observational diversity among SN la and this 
is of prime interest in their application as distance indica- 
tors. Subtle but unequivocal differences between events with 
large wavelength coverage and dense temporal sampling are 
evident in light curve shapes, colour evolutions, luminosities, 
evolution of spectral lines and expansion velocities derived 
from the line shifts. The search for accurate correlations be- 
tween photometric and spectroscopic properties could im- 
prove the luminosity calibration and help to shed light on 
the explosion mecha nisms and progenitor system. Recently, 
iBenetti et all (|2005l ) have identified three classes of SNe la 
based on their spe ctroscopic features. A similar classification 
is also reported in [Branch et al.l (j2006) . Interestingly, two of 
these classes are nearly indistinguishable in some of their 
photometric parameters (e.g. Amis, Phillips et al. 1993); 
yet they show clear spectroscopic differences. High velocity 
gradient (HVG) SNe are characterized by a fast decrease 
in their expansion velocity over time as measured from the 
minimum of the Si II (6355 A) absorption line. On the other 
hand in the low velocity gradient (LVG) SNe group, which 
represents the majority of SN la, the evolution of the Si II 
(6355 A) velocity is smooth, the ejecta expansion is slower 
than in HVG SNe. These spectroscopic differences may or 
may not affect the relations used to calibrate the luminosi- 
ties of SNe la. Subtle differences in colour may exist among 
the different classes and lead to significant bias in redden- 
ing estimates. The existence of spectroscopic families sug- 
gests possible differenc es in the progenitor channels an d /or 
explosion mechanisms. [Branch fc van den Berghl l|l993l ) for 
example reported that SNe characterized by high expansion 
velocities, tend to explode in late type galaxies, pointing to 
progenitors arising from a young population. 

In this paper we present the observations of the HVG 
SN 2002dj carried out by the European Supernova Col- 
laboration (ESC). SN 2002dj (a = 13''13'"00^34, 5 = 
-19°31'08".7, J2000) was discovered in NG C 5018 on June 
12.2 UT (lAUC 7918. lHutchings fc L|200^ ) and classified by 
ESC members as a Typ e la event on June 14.15 UT (lAUC 
7919. lRiello et al]|2002l ). The layout of the manuscript is as 
follows. Observations and data reduction are presented in 
Sect. 2. We describe the reddening estimate in Sect. 3. In 
Sect. 4 we analyze the SN 2002dj optical and IR photometry 
comparing its properties with those of its kinematical twin, 
SN 2002bo. We also investigate whether they are representa- 
tive of the HVG SNe group. Sect. 5 contains the determina- 
tion of the SN 2002dj absolute luminosity and a character- 
ization of the properties of its host galaxy. The optical and 
IR spectroscopy is analyzed in Sect. 6 and Sect. 7, respec- 
tively. The expansion velocities are discussed in Sect. 8 and 
models of the early and late spectra in Sect 9. We conclude 
in Sect. 10 by examining possible physical conditions in the 
HGV SNe ejecta that could explain part of their spectro- 
scopic and photometric behavior. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Instrument settings 

SN 2002dj was observed with seven different instruments in 
the optical and two facilities in the near infrared. Below we 
list the instruments and their main characteristics: 



Table 3. Magnitudes for the local IR photometric sequence in the 
field of SN 2002dj (Fig. [T] right panel). The data were obtained 
during throe photometric nights with Sofl. 



id 


J 


H 


Ks 


2 


12.23 ± 0.03 


11.75 ± 0.03 


11.63 ± 0.03 


2 IR 


14.80 ± 0.03 


14.18 ± 0.03 


13.95 ± 0.04 


7 


13.63 ± 0.04 


13.42 ± 0.04 


13.38 ± 0.04 


11 


14.32 ± 0.04 


13.89 ± 0.04 


13.76 ± 0.04 



• 0.9m telescope located at Cerro Tololo Inter- American 
Observatory (CTIO) equipped with a CCD camera 
(2048x2048, pixel size = 0.40 arcsec pixel"^) and UBVRI 
Johnsons, Kron-Cousin standard filters. 

• 1.0m telescope located at CTIO equipped with A Novel 
Double-Imaging CAMera (ANDICAM; Hawaii HgCdTe 
1024x1024, pixel size = 0.137 arcsec pixeP^) and JHKs 
standard filters. 

• European Southern Observatory (ESO) New Technol- 
ogy Telescope (NTT) located at La Silla Observatory and 
equipped with the ESO Muhi Mode Instrument (EMMI) 
in RILD mode (2xMIT/LL CCD 2048x4096, pixel size = 
0.167 arcsec pixel"^) and BVRI standard filters with ESO 
identification numbers 605, 606, 608 and 610, respectively. 

• ESO-NTT equipped with the Son of ISAAC camera 
(Sofl; Hawafi HgCdTe 1024x1024, pixel size = 0.29 arcsec 
pixel"^) and JHKs standard filters. 

• Danish 1.54m telescope located at La Silla Observatory 
equipped with the Danish Faint Object Spectrograph and 
Camera (DFOSC; MAT/EEV CCD 44-82 2048x2048, pixel 
size = 0.39 arcsec pixel"^) and UBVRI standard filters with 
ESO identification numbers 632, 450, 451, 452 and 425, re- 
spectively. 

• ESO-Kueyen Very Large Telescope (VLT) located at 
Paranal Observatory equipped with the FOcal Reducer and 
low dispersion Spectrograph (FORSl; 2xE2V 2048x4096 
pixel size — 0.2 arcsec pixel" ^) and BVRI standard filters 
with ESO identification numbers -1-34, -1-35, 4-36 and +37, 
respectively. 

• Nordic Optical Telescope (NOT) located at Roque de 
los Muchachos Observatory equipped with the Andalucia 
Faint Object Spectrograph and Camera (ALFOSC; E2V 
2048x2048, pixel size = 0.19 arcsec pixel"^) and UBVR 
standard filters with NOT identification numbers 7, 74, 75 
and 76, respectively, and an interference i filter with number 
12. 

• Isaac Newton Group (ING) Jacobus Kapteyn Telescope 
(JKT) located at Roque de los Muchachos Observatory 
equipped with a CCD camera (SITe2 2048 x 2048, pixel size 
= 0.33 arcsec pixeP^) and UBVRI standard filters with ING 
identification numbers 3, 27, 30, 37 and 44, respectively. 

• ING-Isaac Newton Telescope (INT) located at Roque 
de los Muchachos Observatory equipped with the Wide Field 
Camera (WFC; 4 thinned EEV 2048x4096, pixel size = 0.33 
arcsec pixel" ^) and UBVRi(sloan) standard filters with ING 
identification numbers 204, 191, 192, 193 and 215, respec- 
tively. 
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Figure 1. Finding chart for SN 2002dj (left panel) and the local standard stars (right panel) in a R-band exposure obtained on July 25 
2002 with the Danish 1.54m+DFOSC. 



Table 1. Magnitudes of the local photometric standard stars in the field of 
SN 2002dj (Fig. [T] right panel), obtained during seven photometric nights with 
DFOSC, CTIO 0.9m, WFC and FORSl. 



id 


U 


B 


V 


R 


I 


1 




17.91 ± 0.02 


17.09 ± 0.03 


16.61 ± 0.03 


16.18 ± 0.04 


2 


15.79 ± 0.04 


15.08 ± 0.01 


14.09 ± 0.02 


13.53 ± 0.02 


13.02 ± 0.01 


3 




18.11 ± 0.02 


17.33 ± 0.04 


16.88 ± 0.03 


16.43 ± 0.02 


4 




18.15 ± 0.04 


17.51 ± 0.03 


17.13 ± 0.01 


16.76 ± 0.06 


5 


16.44 ± 0.05 


16.45 ± 0.02 


15.84 ± 0.01 


15.48 ± 0.01 


15.11 ± 0.01 


6 




18.59 ± 0.02 


17.52 ± 0.05 


16.90 ± 0.01 


16.31 ± 0.01 


7 


15.05 ± 0.03 


15.08 ± 0.01 


14.62 ± 0.02 


14.33 ± 0.01 


14.04 ± 0.03 


8 


18.01 ± 0.08 


17.60 ± 0.01 


16.57 ± 0.01 


15.96 ± 0.02 


15.38 ± 0.01 


9 




19.25 ± 0.06 


18.68 ± 0.04 


18.32 ± 0.05 


17.85 ± 0.08 


10 


17.31 ± 0.07 


17.28 ± 0.03 


16.62 ± 0.01 


16.22 ± 0.04 


15.85 ± 0.03 


11 


16.87 ± 0.06 


16.76 ± 0.03 


15.98 ± 0.03 


15.50 ± 0.02 


15.03 ± 0.01 



2.2 Optical and infrared photometry 

For the optical photometric observations basic data reduc- 
tion (bias and flat-field correction) was performed using 
standard routines in irafQ. In Table 1 are reported the 
magnitudes of the local photometric sequence identified in 
Fig. [Tl that we re calibrated against Landolt standard stars 
l|Landoltlll99^ ) in the seven photometric nights marked in 
Table 2. SN photometry was performed using the PSF fitting 
technique. Only in the B and V FORSl images obtained on 
March 25, 2003, was it possible to remove the galaxy. For the 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc, under contract to the National 
Science Foundation. 



other instruments the lack of a suitable image without the 
SN has prevented us from using the template subtraction 
technique. We note that the FORSl B and V band frames 
on which the SN is still relatively bright yield magnitudes in 
agreement to the ones derived from the template subtracted 
images. This gives us confidence that the PSF photometry 
performed around maximum, when the SN signal-to-noise 
ratio is high, was not significantly contaminated by the host 
galaxy background. 

The near-IR data reductions were also performed us- 
ing standard IRAF routines. Dark and fiat-field corrections 
were applied to the scientific frames. For each night a sky 
background image was created by taking the median of the 
dithered science frames. For the ANDICAM images this ap- 
proach lead to an imperfect removal of the galaxy light. 
Nevertheless, the good agreement between the SofI and 
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Table 2. S-corrected optical photometry of SN 2002dj. The photometry obtained using the color equation are in Table lA3l 



UT date 


M.J.D. 


Phase" 


U 




B 




V 






R 






I 




Instr. 


1 /OR /9C\f\'?P 


52439.1 


-10.9 


1 ft 07 + n 04 


16.06 


-|- 


0.02 


15.85 i 


0.02 


15.66 


-|- 


0.02 


15.90 


-|- 


0.03 




1 4 /Ofi 79009 




-10.1 




15.64 




0.02 


15.49 i 


0.02 


15.27 


-|- 


0.02 


15.46 


-|- 


0.02 


F-MMT 


97/0fi/9009P 




2.9 




14.36 




0.02 


14.14 ^ 


0.02 


14.10 




0.02 


14.46 


-|- 


0.02 


WFC 


9K /Ofi /9009 


^94^*^ Q 


3.9 




14.41 




0.02 


14.16 i 


0.02 


14.12 




0.02 


14.48 


-|- 


0.02 


JKT 


9Q /Ofi /9009 




5.0 


^AA^ + n (14 


14.46 




0.02 


14.20 i 


0.02 


14.15 




0.03 


14.51 


-|- 


0.03 


JKT 


'^0/0fi/90n9 


52456 n 


6.0 














14.24 




0.02 








JKT 


09 /07/9009 




8.0 




14.66 




0.09 


14.31 i 


0.02 


14.31 




0.02 


14.61 


-|- 


0.02 


JKT 


O'i /07/9009 




8.1 


1 4 ^^fi + n 04 


14.67 




0.02 


14.27 i 


0.02 


14.30 




0.02 


14.66 


-|- 


0.05 


DFOSC 


04/07 /9009 


52460 n 


10.0 


1 K no + n OA 


14.81 




0.02 


14.36 ± 


0.02 


14.46 




0.02 


14.85 


-|- 


0.05 


DFOSC 


07/07/9009P 


^94fi'^ 


13.0 




15.13 




0.02 


14.56 i 


0.02 


14.61 




0.03 


14.94 


-|- 


0.04 


DFOSC 


OQ /07/9009 


^94fi^ 


15.0 




15.39 




0.02 


14.73 i 


0.02 


14.71 




0.02 


14.97 


-|- 


0.02 


DFOSC 


11/07/2002 


52467.0 


17.0 


16.06 ± 0.04 


15.61 


± 


0.02 


14.85 ± 


0.02 


14.76 


± 


0.02 


14.95 


± 


0.02 


DFOSC 


17/07/2002P 


52473.0 


23.0 










15.12 ± 


0.02 


14.86 


± 


0.02 


14.81 


± 


0.02 


CTIO 0.9m 


16/07/2002 


52473.1 


23.1 


16.75 ± 0.04 


16.29 


± 


0.02 


15.13 ± 


0.02 


14.83 


± 


0.02 


14.89 


± 


0.02 


DFOSC 


19/07/2002 


52476.0 


26.0 


16.99 ± 0.04 


16.54 


± 


0.02 


15.40 ± 


0.02 


15.05 


± 


0.02 


14.85 


± 


0.02 


DFOSC 


24/07/2002 


52479.9 


29.9 




16.89 


± 


0.02 


15.58 ± 


0.02 


15.14 


± 


0.02 


14.83 


± 


0.02 


ALFOSC 


25/07/2002 


52481.0 


31.0 


17.44 ± 0.05 


16.91 


± 


0.02 


15.59 ± 


0.02 


15.18 


± 


0.02 


14.88 


± 


0.02 


DFOSC 


29/07/2002 


52484.9 


34.9 


17.59 ± 0.15 
























ALFOSC 


08/08/2002P 


52495.0 


45.0 




17.41 


± 


0.04 


16.25 ± 


0.02 


15.95 


± 


0.03 


15.77 


± 


0.06 


CTIO 0.9m 


31/08/2002 


52518.0 


68.0 


18.43 ± 0.09 








16.83 ± 


0.02 














DFOSC 


02/09/2002 


52520.0 


70.0 










16.87 ± 


0.02 














DFOSC 


06/09/2002 


52524.0 


74.0 




17.75 


± 


0.19 






17.06 


± 


0.06 


17.01 


± 


0.48 


DFOSC 


25/03/2003P 


52724.3 


274.3 




21.00 


± 


0.04 


20.71 ± 


0.03 


21.39 


± 


0.06 


20.78 


± 


0.09 


FORSl 


22/05/2004P 


53147.2 


697.2 




>24.7'' 


>23.4'' 














FORSl 



° Counted since the time of the B maximum brightness M.J.D. =52450 ± 0.7 

^ Upper limit 

^ Photometric night 



ANDICAM photometry gives us confidence that the small 
galaxy residuals did not significantly bias the ANDICAM 
magnitudes. The sky background image was then subtracted 
from each single frame. For SofI, an illumination correction 
was also applied to all images. Four stars close to the SN 
position were c alibrated in th e JHK s bands using stan- 
dard stars from iPersson et al.l (|l998l ) during three photo- 
metric nights marked in Table |4] The resulting magnitudes 
are listed in Table |3] They agree well with those reported 
in the 2MASS catalog with Jseq - J2mass=0.01 ± 0.01, 
Hseq - -H'2AfAss=0.00 ± 0.01 and A's^eq - Ksimass = 
0.01 ± 0.01. The final SN calibration in all bands was per- 
formed using the S-correction technique (see Appendix for 
details). Thanks to the very good spectroscopic coverage 
we could compute the S-terms using only the spectra of 
SN 2002dj for the B, V, R, I bands, while for the U, J, 
H, Ks, the SN 200 2d,i spectra were com plemented with 
those of SN 2005cf (|Garavini et al.l |2007| ) and SN 2002bo 
ijBenetti et al.ll2004l ). 



1 1994 ) and checked against the photometry. When discrepan- 
cies occurred, the flux of the spectrum was scaled to match 
the photometry. In nights with good observing conditions 
the agreement with the photometry was found to be within 
10%. 

The near-IR spectra were reduced u sing standard pro - 
cedures in the FIGARO 4 environment (jShortridgel 119951 ). 
Wavelength calibration was carried out using Xe arc spec- 
tra and the accuracy of the solution was checked using OH 
sky lines. The flux calibration was secured with respect to 
near-IR solar-analog standard stars observed close in time 
and airmass to the S N to minimize variations in th e at- 
mospheric absorptions (jMaiolino. Rickc fc Rickc"l996!'l . The 
SN spectrum was divided by the standard stars spectrum to 
remove the strong telluric IR features. The resulting spec- 
trum was then multiplied by the solar spectrum to eliminate 
the intrinsic features introduced by the solar-type star. Like 
in the optical, the flux calibration was checked against the 
photometry and, if necessary, scaled to match the latter. 



2.3 Optical and infrared spectroscopy 

Optical spectra were reduced using the ctioslit package in 
IRAF. Optimal extraction was usually obtained by weight- 
ing the signal according to the intensity profile along the 
slit. Sky subtraction was carried out by fitting a low or- 
der polynomial on either side of the SN spectrum and the 
wavelength solution was determined from arc spectra. The 
wavelength calibration was checked against bright night-sky 
emission lines. Flux calibration was p erformed by means 
of spectrophotometric standard stars (|Hamuv et al.l 1 19921 . 



3 INTERSTELLAR EXTINCTION 

Presence of dust in NGC 5018 has be en claimed insev- 
eral studies. The map produced by iFort et all (|l986l ) 
sh ows a dust lane e r nbedd ed in the galaxy. Subsequently , 
iCarollo fc Danzieer I |l994 ) and iGoudfrooii et all (|l994ah 
detected patchy dust in the galaxy core. None of those 
maps shows evidence of dust at the location of SN 2002dj. 
Our deep images taken w- ith the VLT add details to the 
ICarollo fc Danziger I l|l994h colour maps, but again do not 
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Table 4. S-corrected IR photometry of SN 2002dj. The photometry obtained using 
the color equation are in Table lA4l 



UT date 


M.J.D. 


Phase" 




J 






H 




Ks 




Instr. 






-11.0 


15.63 


-|- 


0.06 


15.64 


-|- 


0.09 


15.65 


± 


0.13 




1 79009 


'^9A40 


-10.0 


15.27 


-|- 


0.04 


15.31 


-|- 


0.06 


15.25 


± 


0.09 


ANDTP AM 


1 ^ /Ofi/9009 


52441.1 


-8.9 


14.97 


-|- 


0.03 


15.13 


-|- 


0.03 


15.11 


± 


0.03 


SoFi 




52444.1 


-5.9 


14.63 


-|- 


0.04 


14.84 


-|- 


0.03 


14.71 


± 


0.03 


SoFi 


2n/nfi/2nn2 


^244fi 


-4.0 


14.58 


-|- 


0.11 


14.80 


-|- 


0.07 


14.66 


± 


0.08 




24/nfi/2nn2 

irt / yjyj / ^uvj^ 


52450 


-0.0 


14.60 


-|- 


0.03 


14.84 


-|- 


0.06 


14.52 


± 


0.10 




9il/0fi/9009P 




-0.0 


14.62 


-|- 


0.03 


14.83 


-|- 


0.03 


14.54 


± 


0.03 


SoFi 


27/06/2002 


52453.1 


3.1 


14.86 


± 


0.05 


14.90 


± 


0.04 


14.68 


± 


0.09 


ANDICAM 


30/06/2002 


52456.0 


6.0 


15.30 


± 


0.13 


15.09 


± 


0.10 


14.88 


± 


0.20 


ANDICAM 


08/07/2002 


52464.0 


14.0 


16.30 


± 


0.09 


15.05 


± 


0.05 


15.04 


± 


0.09 


ANDICAM 


11/07/2002 


52467.0 


17.0 


16.26 


± 


0.05 


14.97 


± 


0.03 


14.90 


± 


0.06 


ANDICAM 


11/07/2002P 


52467.0 


17.0 


16.20 


± 


0.03 


14.95 


± 


0.03 


14.83 


± 


0.03 


SoFi 


14/07/2002 


52470.0 


20.0 


16.08 


± 


0.13 


14.81 


± 


0.05 


14.75 


± 


0.08 


ANDICAM 


17/07/2002 


52473.0 


23.0 


15.86 


± 


0.04 


14.75 


± 


0.03 


14.65 


± 


0.05 


ANDICAM 


25/07/2002 


52481.0 


31.0 


15.57 


± 


0.07 


14.79 


± 


0.05 


14.81 


± 


0.09 


ANDICAM 


28/07/2002 


52484.0 


34.0 


15.63 


± 


0.05 


15.06 


± 


0.04 


15.12 


± 


0.07 


ANDICAM 


31/07/2002 


52487.0 


37.0 


15.89 


± 


0.03 


15.23 


± 


0.07 


15.43 


± 


0.29 


ANDICAM 


07/08/2002 


52494.0 


44.0 


16.58 


± 


0.14 


15.61 


± 


0.09 


15.65 


± 


0.35 


ANDICAM 


10/08/2002 


52497.0 


47.0 


16.68 


± 


0.10 


15.73 


± 


0.06 


15.93 


± 


0.16 


ANDICAM 


13/08/2002 


52500.0 


50.0 


16.89 


± 


0.12 


15.85 


± 


0.06 


16.04 


± 


0.15 


ANDICAM 


30/08/2002 


52517.0 


67.0 


18.20 


± 


0.38 


16.66 


± 


0.03 


16.84 


± 


0.03 


SoFi 



" Counted since the time of the B maximum brightness M.J.D. =52450 ± 0.7 
P Photometric night 



show the presence of dust at, or close to, the SN position. 
SN spectra can be used to measure the equivalent width 
{EW) of the Na I D lines as a proxy to quantify absorp- 
tion. The Galactic colour excess we obtain using; one o f 
the two linear relations reported in iTuratto et al.l (|2003l ) 
{E{B -V) = 0.16 X EW yielding E{B - V)GaiacUc = 
0.11 ±0.03 from EW = 0.72 ± 0.19 A) is in very good agree- 
ment with that reported by ISchlegel. Finkbeiner fc Davis! 
l|l998t) {E{B - V)GaiaMc = 0.096) for the SN 2002dj line 
of sight. No absorption line of interstellar Na I D has been 
detected from the host galaxy. Since SN 2002dj was observed 
both in the optical and in the near IR, we can apply several 
methods to estimate its extinction from the observed colour. 
The results are reported in Table [5] We note that all meth- 
ods involving B — V point to an E{B — V) greater than that 
due to Milky Way absorption alone, indicating a small but 
not negligible dust extinction towards SN 2002dj in the host 
galaxy, while the methods based on V — IR colours provide a 
negligible colour excess for the host galaxy. A similar appar- 
ent inconsistency occurred in the case of SN 2002bo, where 
the reddening derived from B — V colour seems to be larger 
than that derived throu gh detailed synthetic spectra mod- 
eling IstehieeEaLlliooi) and, as in the case of SN 2002dj, 
it is significantly larger than that derived from the V — IR 
colours. The peculiar colour s of SN 2002bo were already 
noted bv lBenetti et al.l (|2004l ') and lKrisciunas et al.l (l2004bl ). 

With the aim to check whether HVG SNe could in gen- 
eral have intrinsically different colours which could bias their 
reddening estimate, we examined the other HVG SNe re- 
ported in the literature for which multi-epoch spectra show 
a persistent high expansion velocity (see Table |8]0. Unfor- 

2 Although SN 2002er was classified as HVG SNe llBenetti et al.l 



tunately, all objects except SN 2002bf and SN 2004dt, show 
signs of absorption in their host galaxies (i.e. Na I D lines). 
The latter two SNe are also not useful for this purpose be- 
cause the photometric follow-up for SN 2002bf started ~ 
at +6 days and definitive photometry is not yet available 
for SN 2004dt. Hence, for aU objects in Table |8] the intrin- 
sic colour is rather uncertain. To investigate the possible 
peculiar colours for HVG SNe we compared the the scatter 
between the colour excess derived using the B — V and V — I 
colours among LVG SNe with the systematic difference we 
found for our sample of HVG SNe. 

To determine the differe nce in extinction me asure for 
LVG SNe we used Table 2 of iPhillips erall (|l999l ) and ob- 
tained E{B - V)b-v-E{B - V)v-i = -0.01 ± 0.08. This 
value was derived by selecting SNe with 0.95 < Amis < 1-60 
in order to avoid co lour pecul iarity not accounted b y the 
iPhiUips et aU l|l999t ) relation l|Garnavich et al.1 120041 '). For 
HVG SNe in addition to SN 2002bo and SN 2002di, only in 
the case of SN 1997bp, B FT?/ photometry llAhaviUa et al.1 
2004; Jha et al. 2006i ) was suitable to compute the colours 
at maximum. For these three objects we find E{B — V)b-v- 
E{B - V)v-i = 0.17 ± 0.03. The systematic difference in 
reddening depending on the reference color in HVG SNe is 
hence nearly 2a the scatter observed in LVG SNe hinting 
at intrinsic colour differences. Anyway, the colour evolution 
among HVG SNe shows significant variations (see Fig. [6} 
making it difficult to establish an overall difference with LVG 
SNe. 



|2005^ it was not included in the table because of its lower ex- 
pansion velocities. This SN could indeed represents a transi tion 
object between HVG SNe and LVG SNe llTanaka et al.ll2008h . 
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Figure 2. UBVRI light curves of SN 2002dj. Our data are shown with filled circles using di fferent colours for diff erent instruments. 
For comparison the light curves of other SNe as well as the te mplate for BVI Ami5 = 1.11 by [Hamuy et al. I lll996h (d a,shed lines) are 
displayed. HVG SNe a re indicated b y filled symbols : SN 2002bo llBenetti et al.ll2004lKrisciunas et al.ll2004d )~SN 1997bp (|AItavilla et al.l 
|2004 |jha et al.ll2006l). SN 2002bf jLeonard et al.l 120051 SN 2006X llWane et aLlboOSah , while LVG SNe are drawn with open sym- 
bols: SN 1998bu lljha et al.. Suntzeff et al.lll999l ). SN 2001el l|Krisciunas et al.ll2003h . SN 2005cf l|Pastorello et al.ll2007^ and SN 2003du 
llStanishev et al.ll2007t) . To facilitate the comparison the light curves of the different SNe were normalised to match at maximum (between 
—5 to +5 days). 



4 OPTICAL AND INFRARED PHOTOMETRY 



4.1 Light curves 



SN 2002dj was observed in UBVRI from -11 to +697 days 
and from —11 to +68 days in JHKs. The corresponding 
light curves are shown in Figs. [2l[3l|4l[5l while colour curves 
are displayed in Figs. [6] and [7] 



In Fig.Elthe BVI and, to a lesser degree i?, HVG SNe (filled 
symbols; SN 2002dj, SN 2002bo and SN 2006X) show a faster 
rise to maximum if compared to LVG SNe (empty symbols; 
represented by SN 2001el, SN 2003du and SN 2005cf). On 
the other hand we note that SN 2006X and SN 1997bp peak 
later in VRI (~+2 days) than the LVG average popula- 
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Table 5. SN 2002dj total (Galactic 
from different methods. 



host galaxy) reddening 



Method 



E{B - V) Reference 



Bjnax Vmax 

^max ^max 

(B — V)c MAGIC 

{B - V)jha 

(B - V)ta^l 

(B — V)wang 

{B-V)v-H 

{B-V)v-K 



0.22 ± 0.06 Phillips et al. (1999) 
0.24 ± 0.06 Altav illa et al. (2004) 
0.10 ± 0.06 Phil lips et al. (1999) 



Philh£s_et_aL 

0.32 ± 0.06 Wang et al. (2003) 

0.29 ± 0.06 Jha et al. (2006) 

0.27 ±0.14 Phillips et al. (1999') 

0.21 ± 0.04 Wa ng et al. (2005) 

0.05 ± 0.04 Krisciunas et al. (2004a) 

0.10 ± 0.06 Krisciunas et al. (2004a) 
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Figure 3. Early phase UBVRI light curves of SN 2002dj. The 
symbols, SNe and vertical shifts are the same as reported in Fig. (2] 



tion (see histograms in IContardo et al] (|2000l )). Therefore 
the fast rise in VRI of the HVG objects could be induced 
by a peculiar behaviour of the B hght curve that histori- 
cally is used to set the phase of the SN la. Comparing the 
light curves of SN 2006X (SN 1997bp does not have ob- 
servations close to maximum) with those of SN 2005cf and 
leaving the B maximum time as a free parameter in the 
fit, we found a phase shift of —0.9 days for the V and R 
bands and —0.5 for the /. The latter results reduce the pos- 
sibility that the early time behaviour of HVG SNe could 
simply be a phase mis-calibration. This is also confirmed by 
the IR observation of SN 2002dj, SN 2002bo and SN 2006X 
which show that those SNe reac h the peak at the same epoch 
as the iKrisciunas et al.l l|2004al ) templates (Fig. |3| . For the 
U band the situation is more complex. The data are more 
scattered and systematic errors due to standard filter mis- 
matching could play an important role. Therefore a photo- 
metrically more homogeneous data set is necessary to con- 
firm the trends observed in the other filters. 

Another feature visible in Fig. [2] is that the HVG 



SN 1997bp, SN 2002b£] and SN 2006X show more pro- 
nounced infiections around 4-25 days in V and R. In the 
case of SN 2002bf the large photometric errors associated 
with the corresponding points demand caution, while for 
SN 2006X, at least part of this effect could be due to its 
high extinction shifting the effective wavelength of the pass - 
band to redder l|Ehas-Rosa et al.ll2006l : IWang et aLlbOOSal ). 
The absorption suffered by SN 1997bp is low and the pho- 
tometric errors small, hence this feature must be intrinsic. 
Finally, we note that SN 1997bp, SN 2002bf and SN 2006X 
are the highest velocity SNe la in our sample, with veloci- 
ties about 5000 km s"'^ higher than those of SN 2002dj and 
SN 2002bo (Fig. [21]). This suggest that photometric features 
could be correlated with the kinematics of the ejecta. 

Since, as previously mentioned, at early epochs the hght 
curves of HVG SNe are different from those of LVG SNe, we 
decided to constrain the epoch of the SN 2002dj B maxi- 
mum brightness, hereafter T{Bmax), using only the B hght 
curve of SN 2002bo. In fact, it is not viable to determine 
the maximum brightness directly from the SN 2002dj hght 
curve because it happens to fall in a gap of 12 days with 
no photometric observations. Having constrained T{Bmax), 
(M.J.D. =52450.0 ± 0.7) we used the light curves of SNe 
with similar behavior around maximum to estimate the 
SN 2002dj peak magnitude in UBVRI (Table [gj. The val- 
ues o f the hght curve parameters Amis and stretch fac- 
tor s l|Perlmutter et anil997a ) are reported in Table [7] The 
stretch facto r was obtain e d by fitting the SN 2002dj B hght 
curve to the iLeibundgutI l| 19881 ) template complemented by 
iGoldhaber et all l|200ll )7 We note that th e values of Amis 
and s, satis fy the relations reported by iPerlmutter et all 
l|l997bl 'l and lAhavilla et al.l |200J) reasonably weU. 

As in the optical, in the near infrared bands both 
SN 2002dj and SN 2002bo show a faster rise to m a jcimum 
brightness than SN 2001el and the IKrisciunas et all (|2004ar ) 
"mid-range" decliner IR templates. SN 2006X follows the 
templates closely, although very early time observations are 
missing. 

At epochs >-|-35 days in 7, R and V HVG SNe are 
perhaps slightly brighter, but the points show large spread, 
covering possible subtle systematic differences. In B HVG 
SNe are clearly brighter than LVG SNe. Interestingly, as 
in the case of the inflections occurring around -1-25 days, 
the brightest objects among H VG SNe are t hose displaying 
the highest expansion velocities. IWang et al] (2008a) aheady 
pointed out that the decline rates 7 between -1-40 and -f-lOO 
days (Pskovskii 1984) of the HVG SN 1984A, SN 2002bo 
and SN 2006X are smaller than for "normal" SNe la (i.e. 7 
= 1.40 ± 0.10 mag 100 d"^), being 1.14 ± 0.06, 1.17 ± 0.10 
and 0.92 ± 0.05, respectively. Fitting the three B observa- 
tions of SN 2002bf at phase > -1-40.0 days we obtained 7 
= 0.90 ± 0.60. The large error is due to the noisy photom- 
etry. For SN 2002dj, considering the only two B observa- 
tions in the phase interval, we derived 7 =1.17 ± 0.13. Also 
for SN 1997bp we have only two points from which yield 7 
=0.83 ± 0.18. Therefore, although with limited statis tical 
strength, we confirm the finding of IWang et al.l (|2008al ). In 
U the fastest expanding SN 1997bp and SN 2006X seem to 



^ The light curves of SN 2002bf have been K-corrected to the 
host galaxy rest frame velocity of SN 2002dj 
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Figure 4. JHK light curves of SN 2002dj, different colours 
refer to different instr uments. For comparis on the light 
curves of SN 2002 bo llKrisciunas et al.l l2004bh. SN 2006X 
dWang et all l2008al). SN 1998b u l lHernandez et al.l I2OO0I). 
SN 2 001el l|Krisciunas et aI.|[20o3 ). SN 2005cf llPastorello et al, 
l2007t ) and the mid range decliner templates from lKrisciunas et al. 



(|2004a|) (solid line) are displayed. The light curves of different SNe 
were vertically shifted in order to match them at maximum. 



follow the B band trend, while SN 2002dj and SN 2002bo 
are closer to the evolution of SN 2005cf and SN 1998bu. 

At late phases SN 2002dj, although an HVG SN resem- 
bles the behavior of LVG SNe. At the latest epoch for which 
we obtained measurable data (+274 days), SN 2002dj ap- 
pears in fact in all bands close t o the light curve of SN 2001el 
ijStritzinger fc Sollermanll2007l . see Fig.^J. SN 1997bp is in- 
stead ~0.5 mag and ~0.3 brighter in B and V , respectively. 
Its late time decline rate is similar to the other SNe, which 
argues against a possible contribution of a light echo. With 
the aim to further investigate (or exclude) the presence of 
a light echo, B and V imaging of SN 2002dj was performed 
with FORSl on April 22, 2004, i.e. 697.2 days after B max- 
imum light. Nothing is visible at the position of the SN and 
through re-detection of artificial stars, we estimated 3cr up- 
per limits of B ~24.7 and 1/-23.4. 



4.2 Colour curves 

The de-reddened colour curves of SN 2002dj are compared 
in Fig. [6] with those of SN 1997bp, SN 2002bo, SN 2002bf, 
SN 2001el, SN 2003du and SN 2005cf. For HVG SNe the red- 
dening was removed using the values of E{B — V) reported 
in Table |8] and assuming Rv=S.l. Because of their large red- 
dening, SN 2002bo and SN 2006X are not considered in the 
HVG vs. LVG comparison. 

In the {B — V)o colour, at least until -1-30 days, HVG 
SNe form a redder sequence than LVG SNe. Afterward the 
two groups mix due to the HVG slower decline rate in the B 
band. Note that the latter poses problems if one wishes to 



Figure 5. Same as Fig. [2l showing the late phase BVRI light 
curves of SN 2002dj. 
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Figure 6. De-reddened {B - y)o, {V - R)o and {V - I)o colour 
curves of SN 2002dj. 



apply the Lira relation (iPhillips et al.lll999l) to estimate the 
reddening of HVG (see also IWang et al.ll2008al ). Therefore, 
both maximum light and tail B — V colours could bias the 
reddening estimation of these objects. In {V — R)o there is 
not a clear separation between HVG and LVG SNe, while 
in (V^ — /)o HVG SNe seem to form a bluer sequence with 
respect to LVG SNe, but the significance of the differences is 
red uced by the large cah bration uncertainties of the I band 
(see iPignata et alll2004l ) . 

The {V - R)o and {V - /)o colours of SN 2002bo are 
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Figure 7. De-reddened (V - J)o, (V - H)o and {V - K)o colour 
curves of SN 2002dj. 

clearly bluer then those of SN 2002dj. This reflects the lack 
of flux in the red part of the SN 2002bo spectrum noticeable 
in the comparison with SN 2002dj (see section 6.3). 

The {V - J)o, {V - H)o and (V - K)o colour s curve s 
of SN 2002dj (Fig. [T]) follow the iKrisciunas et all (|2004al ) 
"mid range" templates reinforcing the idea of negligible host 
galaxy reddening. As in the case of the {V — R)o and {V — I)o 
colours, SN 2002bo is bluer at all epochs. The difference to 
SN 2002dj increases toward redder bands, suggesting that 
the lack of flux starting around 6500 A continues monoton- 
ically towards longer wavelengths. 



Optical and IR observations of SN 2002dj 9 

Table 6. Peak magnitudes of SN 2002dj 



Filter 


m(obs)" 


M(cor)'' 


M(ave)'= 


B 


14.30 ± 0.04 


-19.03 ± 0.23 


-19.34 ± 0.02 


V 


14.15 ± 0.04 


-19.08 ± 0.19 


-19.26 ± 0.02 


R 


14.10 ± 0.05 


-19.06 ± 0.17 


-19.26 ± 0.02 


I 


14.35 ± 0.06 


-18.72 ± 0.14 


-19.00 ± 0.02 


J 


14.56 ± 0.02 


-18.46 ± 0.11 


-18.57 ± 0.03 


H 


14.80 ± 0.03 


-18.18 ± 0.11 


-18.24 ± 0.04 


Ks 


14.51 ± 0.05 


-18.45 ± 0.11 


-18.42 ± 0.04 



° Apparent magnitude 

' Absolute magnitude corrected for reddening 
Average absolute magnitude for SNe la 




20 -10 60 

Days since B maxliriuiii 



5 ABSOLUTE LUMINOSITY 



The Virgo infall model of iKraan-Kortewed (|l986l ) yields a 
distance to relative Virgo of 2.52 for NGC 5018. A ssuming 
a Virgo distance of 15.3 Mpc (jPreedman et al ]|200l') we ob- 
tain 38.55 Mpc (^i = 32.93) for the host gala xy of SN 2002dj. 



The Nearby Galaxies Catalog (ITuIIvI r988l) reports a larger 



distan ce (40.9 Mpc, /i=33.1), while iHilker &: Kissler-Patid 
1 19961 ) using the globular cluster luminosity function de- 
rived a distance modulus ranging between 32.8 and 34.0. 
For a consistent comparison to SN 2002bo we decided 
to assume 38.55 Mpc as the best distance estimate for 
SN 2002dj. With this distance and the reddening esti- 
mate we derived the corresponding absolute magnitudes in 
UBVRIJHKs. The values are reported in Table [S] along 
with th ose obtained from the Amis abs olute magnitude rela- 
ti ons oflPrieto. Rest fc Suntzeff I ( 20061) for BVRI and those 
of lKrisciunas. PhiUips fc Suntzefj (|2004l ') for the JHKs. The 
absolute magnitudes in all the bands are in agreement within 
2a with the values predicted by the previously mentioned re- 
lations, although SN 2002dj appears slightly underluminous. 

Using the UBVRIJHKs observations ofSN2002dj_and 
adding the UV contribution determined bv lSuntzefjil 19961) . 



Figure 8. Comparison between the uvoir light curves of 
SN 2002dj and SN 2002bo. The error bars account only for the 
uncertainties in reddening and photometry. The best fit model of 
SN 2002dj (56Ni=O.45M0) is also shown. 



we constructed a uvoir light curve that is compared in Fig. [8] 
with that of SN 2002bo and a model generated us ing a 



Monte Carlo light curve code (|Cappellaro et al 



curve of SN 2002bo published bv lBenetti et al 



I1997D The 



(|2004l ) has 

been rescaled to the new distance estimate {/i = 31.76 ± 
0.07) reported by Ehas-Rosa et al. (2008). Within the er- 
rors the absolute luminosities of the two SNe are compara- 
ble. The model provides a good flt to the brightness and 
width near the peak, yielding a ^^Ni mass of 0.45 ± 0.04 
Mq for SN 2002dj, but after -f20 days it overestimates the 
observed flux. This is probably due to sudden changes in 
the opacity at advanced stages which our model does not 
properly take into account. 

5.1 An intriguing host galaxy 

NGC 5018 is the dominant giant elliptical of a small group. 
It is peculiar in several respects. First, ISchweizer et al.l 
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Table 7. Main parameters of SN 2002dj and its host galaxy. 



Table 9. Optical and IR spectroscopic observations of SN 2002dj. 



UT date M.J.D. 



Host galaxy 
Galaxy type 
RA (2000) 
Dec (2000) 
Heliocentric 

Recession velocity [km s~^] 
Distance modulus 
EiB - V) 

Date of B max (MJD) 

Offset from the nucleus 

Ami5(B) 

stretch factor in B 



NGC 5018 
E3 

13''13™■01^7 
-19°31'12".8 

2816 ± 1 t 

32.93 ± 0.15 

0.096 it 0.05 

52450.0 ± 0.7 

8". 9 West, 2". 8 South 

1.08 ± 0.05 

0.97 ± 0.02 



epoch" 
(days) 



range 

(A) 



res. 

(A) 



Instr. 



t l lRothberg fc Joseph|[20o3) 



Table 8. Host galaxies of HVG SNe considered in the paper 



SN 


Galaxy 


Type 


Ami5(B) 


E{B - V) 


SN 1981B 


NGC 4536 


SAB 


1.11 




SN 1983G 


NGC 4753 


10 


1.37 




SN 1984A 


NGC 4419 


SB 


1.21 




SN 1989A 


NGC 3687 


SAB 


1.06 




SN 1997bp 


NGC 4680 


Sp 


1.00 


0.18(0.04)" 


SN 2002bf 


CGCG 266-031 


SB(r)b 




0.08(0.04)'' 


SN 2002bo 


NGC 3190 


Sap 


1.17 


0.38(0.10)= 


SN 2002dj 


NGC 5018 


E3 


1.08 


0.096(0.05) 


SN 2004dt 


NGC 799 


SB(s) 


1.21 




SN 2006X 


NGC 4321 


SABbc 


1.31 


1.41(0.04)'* 


" average value between E(B 


-V)b-v 


and E{B - V) 


V-I 



' Prom 
Prom 
From 



Stehle et al.l ||2001 



Wang et al.l ll2008al ') 



classified it as one of tiieir best candidates for a 
recent major merger. Second, although it is morphologi- 
cally classified as gE, its nuclear optical spectrum distin- 
guishes itself by having the weakest measurement of the ab- 
sorption line index Mg2 (0.218) f or its velocity dispe rsion 
among over 400 gEs surveyed by iDavies et al.l l|l987l ) . Its 
UV (lUE) spectrum lacks the prominent UV-upturn short- 
ward of 2000 A, whi ch is typical of old, metal-rich spheroids 
l|Bertola et al1ll993l). Nevertheless , through a spectral index 
study. iLeonardi fc Worthevl l|2000l ) found indications for the 
presence of a relatively young stellar population (~2.8 Gyr) 
in the central region of NGC 5018 with nearly solar metal- 
licity. A similar age (~3Gyr) is propose d by iBuson et al] 
(|2004l ) through a UV spectroscopy study. iKim et al.l (|l988l ) 
detected a H I gas bridge connecting NGC 5018 with 
the nearby spiral NGC 5022, indicating an ongoing flow 
toward the giant elliptical. Furthermore, possible young 
(~ 10^ years) globu l ar clu sters have been claimed by 
iHilker fc Kissler-Patid (1 19961) . Final ly the Ha-|-[N II] maps 
reported in iGoudfrooii et all l|l994bl ) reveal the presence of 
an extended emission distributed as a strongly warped disk 
covering the SN 2002dj position that the authors suggest to 
be associated with star forming regions. 
[Branch fc van den Berghl (|l993l ) reported that SNe la char- 
acterized by high expansion velocities, tend to explode 
in late type galaxies and suggested that those objects 
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3300-9300 
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DPOSC 
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EMMI 
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02/06/18 
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3300-9300 
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DPOSC 


02/06/19 


52444 


1 


-5 


9 


9400-25400 


21 


SofI 


02/06/20 


52446 





-4 





3200-9200 


9 


DPOSC 


02/06/21 


52447 





-3 





3300-9300 


9 


DPOSC 


02/06/24 


52450 











9400-16500 


21 


Son 


02/07/03 


52458 


7 


8 


7 


3400-9300 


9 


DPOSC 


02/07/04 
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10 





3300-10300 


9 


DPOSC 


02/07/07 


52463 





13 





3200-9100 


9 


DPOSC 


02/07/11 


52467 





17 





3300-10300 


9 


DPOSC 


02/07/11 


52467 





17 





9400-25200 


21 


SofI 


02/07/16 


52473 


1 


23 


1 


3300-9100 


9 


DPOSC 


02/07/25 


52481 





31 





3700-10300 


9 


DPOSC 


02/09/02 


52520 





70 





3400-9100 


9 


DPOSC 


03/02/01 


52672 


3 
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could have younger progenitors than "normal" SNe la. The 
HVG SNe analyzed in this paper confirm this (Table |8]). 
iHamuv et al.l ()2000l ) pointed out that slow declining events 
preferentially occur in late type galaxies, while fast decliners 
oc cur in all type of gal axies. A similar result was obtained 
bv lSuUivan et al.l (|2006l ) for a sample of high-z SNe la. These 
observations seem at odds with the fact that SN 2002dj was 
hosted by an elliptical galaxy. However, given the peculiar- 
ities of NGC5018 it is possible that SN 2002dj could be 
associate with a relatively young stellar population. 



6 OPTICAL SPECTRA 

Spectra of SN 2002dj have been obtained at 20 epochs span- 
ning phases between —10.9 and +274 days (Table^, follow- 
ing the rapid evolution during the early epochs in detail and 
sampling the late phases more sparsely. The early discovery 
of SN 2002dj allows us to compare its spectral evolution 
(Fig. H) with that of SN 2002bo and other SNe, especially 
at early epochs. 



6.1 Premaximum phase 

In Fig.[T0]we show the first spectrum of SN 2002dj obtained 
about 11 days before B maximum together with coeval spec- 
tra of SN 2002bo, SN 2003du, SN 2005cf and SN 1994D. 
Already at a first glance it is evident how the minima of 
the Ca II H&K, S II and Si II absorption features are more 
blueshifted in the spectra of SN 2002dj and SN 2002bo than 
in LVG SN 2003du, SN 2005cf and SN 1994D. The Si II 
lines are also deeper in HVG SNe. The blueshift is due 
to the higher expansion velocities of the HVG objects (cf. 
Fig. [2lT). The blend ing of high velocity detached features 
(|Branch et al.l [20041 ) with the photospheric line component 
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Figure 9. Optical spectral evolution of SN 2002dj. With the ex- 
ception of the last three spectra, which for presentation purposes 
have been multiplied by a constant, the other spectra of the se- 
quence have been only shifted vertically. The ffi symbol shows 
the position of the main telluric features. The spectra are labeled 
with the epoch in days past B maximum. 
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Figure 10. Spectrum of SN 20 02dj taken at -10. 9 days. The 
coeval spect r a of SN 2002bo l lBen etti ct al.' '2004|) SN 1994D 
iPatat et all Il996l'). SN 2003d u jstanishcv et _alj |2007| ) and 
SN 2005cf l lGaravini et al.ll200'^ are shown for comparison. The 
spectra have been corrected for reddening and redshift. 



could also play an important role in blueshifting the ab- 
sorption minimum of Ca II H&K and Si II A6355 A in 
HVG SNe as it is for the Ca II IR triplet (Mazzali et al. 
2005, Tanaka et al. 2008). For example, the minimum of 
Si II A6355 A in SN 2002dj and SN 2002bo is at the 
same position o f the high velocity component identified by 
iGaravini et~al] (|200if ) in the spectrum of SN 2005cf. The 
feature around 4400 A marked in Fig. \W\ is attributed to 
Si III AA4553,4568 A. The strength of this line correlates 
with temperature and is clearly detected in LVG SNe. The 
spectrum of SN 2002dj only has a hint Si 111, and it is com- 
pletely absent in SN 2002bo. As shown in Fig. 1151 the low 
temperature of HVG SNe at early epochs is confirmed by 
their high i?(Si II). If well traced by i?(Si II) the tempera- 
ture should rapidly rise in HVG SNe and quickly reach that 
of LVG SNe. Indeed at -4 days (Fig. the Si III is vis- 
ible with similar strength in all SNe but SN 2002bo. The 
latter SN also lacks the weak absorption at ~5550 A also 
attributed to Si 111 (A5740 A) that is present in all other 
SNe. It is interesting to note how the minimum of the Si III 
line is at nearly the same wavelength for all SNe, indicat- 
ing that the line is formed at roughly the same position in 
velocity space. The Ca II H&K doublet (blended with Si 11 
A3858 A) appears to be double-dipped in LVG SN 1994D 
and SN 200 3du, while it is a sing le feature in SN 2002dj and 
SN 2002bo. iHatano et al] (|2000l i suggested that this differ- 
ence is related to higher expansion velocities that lead to 
a greater degree of line blending. The Ca II line blending 
could also be enhanced by the presence of h igh velocity fea- 
tures. On the other hand lLentz et all l|2000l ) argued that the 
double-dip disappears when Si II is weaker than Ca II. Given 
that in all fast expanding SNe the Si 11 lines are stronger 
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Figure 11. Same as Fig. 1101 but for a phase of —4 days. The 
bibliographic sources for the spectra of SN 2002bo, SN 1994D, 
SN 2003du and SN 2005cf are the same as those in Fig. [TO] 



than in the slow expanding SNe, line blending seems a most 
consistent explanation. 
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6.2 Postmaximum phase 

A week after maximum the Si II A6355 A line in the 
HVG SN 2002dj and SN 1997bp remains blueshifted with 
respect to the LVG SN 1994D, SN 2003du and SN 2005cf 
(Fig. [T2j. The Si II A6355 A feature is still deeper in 
SN 2002dj and SN 1997bp than in LVG SNe, while, con- 
trary to the pre-maximum spectra, the S II line A5640 A is 
weaker. SN 2002dj and SN 1997bp also show less substruc- 
ture in the region between 4700 and 5100 A with re- 
spect to SN 1994D, SN 2003du and SN 2005cf. Finally we 
note that the absorption minimum of the feature around 
5700 A at this epoch attributed to a blend of Na I D and 
Si II A5972 A (|Garavini et all |2007| ) has about the same 
position, but in SN 2002dj and especially in SN 1997bp 
it is weaker than in the other SNe. By one month after 
maximum the spectra are dominated by iron peak elements 
(Fig. I13|l . The absorption minimum of the lines are nearly 
at the same position in all selected SNe. The early fast ex- 
pansion of SN 2002dj, SN 2002bo and SN 1997bp is now 
only reflected in the broadness of the Ca II IR triplet. The 
Na I D line is still weaker than in SN 2003du and SN 2005cf 
and is now slightly rcdshifted. The line at ~5300 A iden- 
tified by Branch ot al. (2005i) as due to Cr II is barely vis- 
ible in the faster expanding SNe while it remains visible in 
the slower ones. Finally, around two months after maximum 
brightness the features due to Fe II are slightly redshifted 
or at the same position in SN 2002dj and SN 2002bo with 
respect to SN 2003du and SN 2005cf (Fig. [TH). Also the 
Na I D line is slightly redshifted and its intensity is compa- 
rable in all SNe. Interestingly, in the region between 4700 
and 5100 A SN 2002dj and SN 2002bo still show less sub- 
structures. At these epochs line blending due to the high 
velocity should not be severe in HVG SNe, so the lack of 
lines points to a temperature or abundance effect. 



6.3 SN 2002dj versus SN 2002bo, a close 
comparison 

We compare the spectra of SN 2002dj with those of 
SN 2002bo at four different epochs (Fig. I16p . In the wave- 
length interval covered by the B and V filters the sim- 
ilarity between the two objects is evident, but around 
6500 A SN 2002bo shows less fiux than SN 2002dj. Since the 
spectra of both SNe were calibrated in different epochs with 
different instruments, we can rule out an instrumental mis- 
calibration. Alternatively the difference could be attributed 
to an incorrect estimate of the reddening of the two objects. 
However, in this case we would expect a gradual deviation 
of the two spectra, but the discrepancy looks more like a 
break possibly due to a lack of opacity in the photosphere of 
SN 2002bo. For the comparison we chose to obtain the best 
fit between 4000 and 6000 A. If instead we scale the spectra 
of SN 2002bo in order to get the best match with the redder 
part of those of SN 2002dj, we end up with a flux excess in 
SN 2002bo at blue wavelength. While this is possible, the 
comparison of the SN 2002bo colours with those of other 
LVG and H VG SNe and t he res ults of the spectral modeling 
reported bv lStehle et al.l (|2005l ) point to a fiux deficit in the 
red part instead of an excess in the blue. 
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Figure 12. Same as Fig. 1101 for a phase of +9 days. The biblio- 
graphic sources for the spectra are the same as those in Fig. 1101 
and Fig. 1111 while the spectrum of SN 1997bp comes from the 
Asiago archive (unpubUshed). 
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Figure 13. Same as Fig. llOl for phase -|-31 days. The bibliographic 
sources for the spectra are the same as those in Fig. 1111 

7 IR SPECTRA 

Although the sample of published IR spectra of SNe la is 
growing rapidly, their number is still small compared with 
those obtained at optical wavelengths. In this context the 
spectral evolution of SN 2002dj in the IR from day -8.9 to 
+ 17 presented in Fig. [l7]constitutes an important contribu- 
tion to the global database. 

Fig. [18] compares the first SN 2002dj IR spectrum 
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Figure 14. Same as Fig. 1101 but for pliase +70 days. The bibli- 
ograpiiic sources for the spectra are the same as those in Fig. 1101 
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Figure 15. Evolution of /?(SiII) for a sample of HVG SNe (blue 
filled symbols) and LVG SNe (empty green symbols). 



with those of SN 2002bo, SN 1994D and SN 1999ee at 
similar epochs. The IR spectrum at this phase is nearly 
featureless with the remarkable exception of an absorp- 
tion at '^lOSOO A. The iden tification of this lin e remains 
open. It was attributed by iMeikle et al.l (|l996l ) to He I 
A10830 A or Mg II A10926 A based on an Local Ther- 
modynam ic Equilibrium (LTE) t reatment. Using a NLTE 
approach, iMazzali fc Lucvl (|l998l ) discussed the conditions 
under which the He line could form, and found that its time 
evolution was at odds with the observed one. They there- 
fore favoured e ither Mg or Si . A M g II identification was 
also claimed bv I Wheeler et all l|l998l ). 

SN 1999ee shows an emission line at ~ 12300 A which 
iRudv et al.l l|2002l ) identified as due to Fe HI in the spectra 
of SN 2000cx. The presence of this feature in those SNe is 
due to their higher photospheric temperature. As expected, 
it is absent in the spectra of both SN 2002dj and SN 2002bo. 

The -5.9 day spectrum of SN 2002dj shows a feature 
around 12100 A (Fig. ^ which was not visible in the 
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Figure 16. Comparison of the reddening and redshift corrected 
spectra of SN 2002dj (solid lines) and SN 2002bo (dotted lines). 
The epochs labeled in the plot are for SN 2002dj; the SN 2002bo 
spectra are coeval within 1.0 day. The SN 2002bo spectra were 
scaled in order to match those of SN 2002dj between 4000 and 
6000 A. 
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Figure 17. Infrared spectral evolution of SN 2002dj. The spectra 
are corrected for reddening and redshift. For graphic exigence the 
spectra have been multiplied by a suitable factor and vertically 
shifted. 
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Figure 18. IR spectrum of SN 2002di taken at -8. 9 days. The 
coeval spectr a of S N 2002bo l lBenetti et al.l |2004|) SN 1994D 
llMeikle et al.l Il996l). SN 1999 ee llHamuv et al.l |2002|) and 
SN 2003du llStanishev et al.ll2007l ') are shown for comparison. The 
spectra have been corrected for reddening and redshift. 



Figure 19. Same as Fig. [18] for phase -6 days for SN 2002dj, 
SN1994D and SN 2003du. The bibhographic sources for the spec- 
tra are the same as those in Fig. 1181 while the spectrum of 
SN 2003cg was published in lElias-Rosa et al.l l l200d) . 



earliest spectrum (cf. Fig. [17}. iMarion et al.l l|2003l ) iden- 
tified this feature as a blend of several Ca II lines between 
12430 A and 12830 A. Like for the Mg II line, the blueshift 
of the Ca II minimum is larger in SN 2002dj than in all other 
SNe. Finally at ~20800 A a weak emission is visible in the 
spectra of S N 2002di. SN 2003du and SN 2003cg attributed 



to 



Si III by I Wheeler eTal 



Following iRudv et al. 



11998). 

1 20021 ) we fitted the two earliest 
IR spectra of SN 2002dj to estimate the black body tem- 
perature. We used only sections of the spectrum free of tel- 
luric features or absorption lines. Assuming E{B — V)— 0.096 
for SN 2002dj we obtain T=11900 K at -8.5 days com- 
pared with T= 12900 K and T=23500 K for SN 1994D and 
SN 1999ee reported in Fig. [18l For the latter two SN e we 
assumed E{B - Vl^O .Oe ± 0.02 ijPatat et al.l Hgg^ ) and 
E{B-V)=0.30 ± 0.04 (|Stritzinger et al.ll2002l )Trespectivelv. 
The high temperature of SN 1999ee (Ami5(B) = 0.96±0.02) 
is consistent with the presence of the Fe III emission feature 
in its spectrum. It is also similar to t he 25000 K temper- 
ature of SN 2000cx l|Rudv et al.ll2002l ') taken at -8.0 days, 
which shows a stronger Fe III emission than SN 1999ee. This 
gives us confidence that the technique provides reliable re- 
sults, when applied to early time spectra. The sequence of 
increasing temperature from SN 2002dj to SN 1999ee con- 
firms the results obtained from _R(Si II) (Fig. I15|l . Indeed, 
for the SN 2002dj at -5.9 days we found a much higher tem- 
perature (T=20500 K). This confirms its fast rise as traced 
by the -R(Si II) and the appearance of a Si III line in the K 
ban d in the spec t rum o f SN 2002dj. 

iRudv et al.l (I2OO2I ') proposed that the black body IR 
spectrum fit can be used to derive also an estimate of ab- 
sorption. We note, however, that the resulting surface 
for temperature and colour excess is nearly fiat providing 



only weak constraints. At -1-17 days we see that several 
lines populate the IR spectrum of SN 2002dj and those of 
other SNe (Fig. I20|) . In particular, a broad P-Cygni fea- 
ture attributed to Fe II (|Marion et al.l 12003) is observed at 
~12350 A. In SN 2002dj the blueshift of its minimum is 
roughly the same as in SN 2003cg and SN 2002er, while 
SN 1 999ee displays higher velocities. Finally, iMarion et al.l 
((2003 ) suggested that the minimum around ~15000 A de- 
fines the transition from partial to complete silicon burn- 
ing. The fact that the velocities inferred by the position of 
this minimum is the same in SN 2002dj, SN 2003cg and 
SN 2002er indicate that these SNe have an iron core of sim- 
ilar size. 



8 EXPANSION VELOCITIES 

In Fig. [21] we show the expansion velocities measured from 
the blueshift of the minima of the Si II (A6355 A), S II 
(A5460 A) and Ca II (H&K). It is evident that at early 
epochs the HVG SN 1984A, SN 1997bp, SN 2002bo and 
SN 2002dj are faster than LVG SNe in all those lines. The 
two groups are separat ed most strongly in Si II v elocity. Nev- 
ertheless, as noted by Ijeffery fc BranchI (I199G| ) strong fea- 
tures like Si II (A6355 A) reliably trace the photospheric ve- 
locity only at early epochs, when the mass above the photo- 
sphere is small. For HVG SNe hke SN 2002dj and SN 2002bo, 
where Si II is particularly strong an d where high velocit y fea- 
tures could play an important role (|Mazzali et a]|[2005l) . this 
is unlikely even at very early times. In SN 20 02bo for ex- 
ample the photospheric evolution computed by lStehle et al.l 
(|2005l ) even at -13 days has a velocity ~2000 km s"'^ lower 
than that measured from the minimum of the Si II line (cf. 
Fig. [21] upper panel). The relatively weak S II (5460 A) 
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Figure 20. Same as Fig. 1181 for pliase +17.0 days. The biblio- 
graphic sources for the spectra are the same as those in Fig. 1181 
and Fig. | 19 | The s pectrum of SN 2002er that was published by 
iKotak et ad ll2005h . 



Figure 21. Velocity evolution of Si II (top panel), S II (mid- 
dle panel) and Ca II (bottom panel) in SN 2002dj, SN 2002bo, 
SN 1997bp, SN 1984A, SN 2002bf, SN 1998bu, SN 2003du, 
SN 2005cf, SN 1994D and SN 2003cg. The solid lines represent 
th e evolution of the ph otospheric velocity of SN 2002bo computed 
bv lStehle et al.l ||2005|). w h ile th e dashed lines show the velocities 
predicted bv lLentz et al model for the case of xlO (top 

line at epoch 0), x3 (middle line) and xl/3 (bottom line) C+O 
solar metallicity. 



seems to be a better tracer of the photospheric velocity. In 
HVG SNe this feature is at higher velocities than in LVG 
SNe, but due to its steep velocity decline it approaches the 
LVG group already before maximum. If we take S II as a 
proxy for the photospheric velocity, its evolution may indi- 
cate that at very early epochs the photosphere in HVG SNe 
is located at larger radii than in LVG SNe, but moving in- 
ward more quickly, it reaches a comparable position already 
before maximum. For Ca II (H&K), some LVG SNe closely 
resemble the velocity evolution of HVG SNe, making this 
feature unsuitable to separate the two groups. 

In the case of Mg II (10926 A), the expansion veloc- 
ity measured from the absorption minimum at -8.5 days is 
-15000 km s"^ for SN 2002dj and SN 2002bo and signifi- 
cantly slower (~11000 km s"^) for SN 1994D. The feature 
appears deeper in SN 1994D than in SN 2002dj (Fig. [22)) . 
We do not observe a change in the line velocity from the 
first to the second epoch in both SNe and infer that even 
earlier than one week before maximum the photosphere is 
below the magnesium layer. Sinc e Mg II is produced in the 
outer layers by oxygen burning (| Wheeler et al.iri998l ). this 
means that only IR spectra observed earlier than this may 
tell us about the boundary between explosive carbon and 
oxygen burning. Nevertheless, we note in Fig. [22] that the 
blue wing of the Mg II line extends up to ~ 19000 km s~^ 
in SN 2002dj and only up to -15000 km s"^ in SN 1994D. 
Hence, the burning front must have propagated further out 
in SN 2002dj than in SN 1994D. 
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Figure 22. Comparison between the Mg II A10926 feature of 

SN 2002dj (solid line) and SN 1994D (dashed line) at 8.5 

days (left panel) and -6.0 days (right). 
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9 SPECTRAL MODELLING 
9.1 Early phase 

We modeled the earliest spectrum of SN 2002dj obtained 
10.9 days before B maximum using the Montecarlo code de- 
scribed in Mazzali & Lucy (1993), Lucy (1999) and Mazzali 
(2000), m odified to in c lude abundance stratification as de- 
scribed in lStehle et alj ('2005). While a good fit to the overall 
spectrum can be obtained for the input parameters reported 
in the caption of Fig. | 23[ it is not possible with the d ensity 
structure of W7 d Nomoto. Thielemann fc Yokoilll984 ) to re- 
produce the high velocity features. One modification that we 
adopted here to overcome this is to increase the mass at high 
velocity. In particular, we assume the presence of 0.04Mq 
of material at velocities above 18000 km s^^. This may be 
caused b y interaction o f the o uter ejecta with circumstellar 
material. iTanaka et al.l l|2008h discuss this and other possi- 
bilities in detail. Let us note that the model suffers only to 
a minor degree from the problem of overfitting the fiux to 
the red of the Si II (A6355 A) line suggesting high line opac- 
ity in all regions. This can not be the case for SN 2002bo 
where the observed spectrum is noticeably below the model 
l|Stehle et al.ll2005D . 

Finally, we note that the spectral modelling seems to 
exclude the presence of C II in the earliest spectrum of 
SN 2002dj. The presence of this high velocit y feature was 
found in the -9 days spectrum of SN 1998aq (|Branch et al.l 
l2003l'l and in the sp ectra of SN 2006gz from -14 to -10 days 
l|Hicken et al]|2007l ). High velocity C II (A6578 A ) in the -14 
days s pectrum of SN 1990N was also cl aimed bv, Fi sher et al.l 
([l993), although iMazzali et al] l|200ll ) showed that a high- 
velocity component of Si II (A6355 A) was another viable 
identification. Moreover, the presence of small quantities of 
carbon has been show n to improve t he fit of the very early 
spectra of S N 1994D dBranch et aLir2005i ). SN 2001el and 
SN 2003du (|Tanaka et all 120081 '). The latter SN la all are 
LVG SNe, therefore the absence of C II in the -10.9 days 
spectrum of SN 2002d,j and in the -12.9 days spectrum of 
SN 2002bo fetehle et a l 2005) could indicate that the car- 
bon burning has p enetrated further out in these explosions 
l|Stehle et al.ll2005l ). 



9.2 Nebular phase 

In order to estimate the errors on the ^^Ni mass esti- 
mate obtained by modelling nebular spectra, we gener- 
ated a set of synthetic spectra computed using two dif- 
ferent codes and two different assumptions of the ^^Ni 
mass distribution. The six synthetic spectra together with 
the two observed spectra of SN 2002dj are presented in 
Fig. 1241 Both codes assume line-formation in non-local ther- 
modynamic equilibrium ( n on-LT E) conditions. The code 
described in iBowers et all l| 19971 ) assumes a uniform den- 
sity, homologously expanding sphere containing iron, cobalt 
and sulphur, with relative abundances modulated by the 
radioactive decay. In computing the SED only singly and 
doubly- ionized species are considered because they are pre- 
dicted to dom i nate the spectrum. The code described in 
iMazzah et all l|200ll ) can use either a homologously ex- 
panding nebula of finite exten t , unif orm density and com- 
position, as in iBowers et al.l (|l997l ). or stratification in 




1 [A] 

Figure 23. Comparison between the first spectrum of SN 2002dj 
taken at -10.9 days and a synthetic spectrum. The the input pa- 
rameters for the latter are E{B — y)=0.096, Vpji=11600 km s~^, 
log(L)=42.33. 



density and abundance ([Mazzali et al.l[2007h . In this case, 
the density profile is taken from th e explosion model W7 
iNomoto. Thielemann fc Yokoi|[l983 '). 7-rays and positrons 
are emitted at various depths according to the distribution 
of ^^Ni and their propagation and deposition is followed 
using a Monte Carlo sc heme similar to that discussed by 
ICappellaro et af] (|l997[ l for their light-curve models. Gas 
heating and cooling as well as line emissivity are computed 
in non-LTE in each radial shell, and the line profiles are 
computed assuming each shell to contribute to a truncated 
parabola, with an inner truncation point corresponding in 
velocity to the inner boundary of the shell considered. The 
sum of these truncated parabolas is the emerging spectrum. 
Constant 7-ray (k^, — 0.027 cm'^ g~^) and positron opacities 
(Kg+ = 7cm'^g^^) are assumed in both cases. 

The synthetic spectra yield slightly different results for 
the ^®Ni mass, reflecting t he different spectral calibrations 
and code assumptions. The IBowers et all (| 19971 ) model gives 
a '^''Ni mass of 0.59 Mq and 0.55 Mq for the -|-222d and 
-|-274d spectra, resp ectively, while the one-zone model from 
iMazzali et"aLl (120011) yie lds 0.55 Mq for both epochs. With 
the Mazzali et al. f (|2007l ) stratified model we obtain 0.45 Mq 
for both spectra. Results obtained with the same code on the 
two different epochs are consistent, indicating the robust- 
ness of the fit. On the other hand, both one-zone models 
yield a '^'^Ni mass slightly larger than the stratified model. 
iRuiz-Lapuente et al.l l|l995h found similar results with their 
stratified nebular model as o pposed to the one-zone model 
(|Ruiz-Lapuente fc Lucvl[l992l ) . The main difference between 
the models is that in the one-zone codes the ^®Ni density is 
uniform, while in the stratified code it is higher near the cen- 
tre, which leads to a higher 7-ray deposition at late times 
and reduces the ^^Ni mass necessary to produce a given lu- 
minosity. The ®^Ni mass obtained with the stratified model 
is in excellent agreement with the one obtained from the 
uvoir light curve, suggesting that this assumption provides 
a good description of the physical conditions in the nebula. 
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Figure 24. Nebular spectra of SN 2002dj taken at +222 days (top 
panel) and +274 (bottom panel) compared with various synthetic 
spectra. 

10 DISCUSSION AND CONCLUSION 

We have presented optical and IR photometric and spec- 
troscopic observations of the HVG Type la SN 2002dj. The 
ejecta kinematics of this object are very similar to that of 
SN 2002bo. The direct comparison of the line shapes and 
strengths of the two objects reveals an even more impres- 
sive similarity. The one obvious distinction is a peculiar flux 
deficit redward of 6500 A in SN 2002bo. 

The light curves are nearly identical with the notable 
exception of the B band at epochs later than +40 days. 
However, the V — IR colours are remarkably different. 
SN 2002bo is much bluer confirming the lack of fiux redward 
of ~6500 A. Finally, also the bolometric absolute luminos- 
ity around maximum of the two objects turns out to be very 
similar. The IR contribution is in in fact very small at those 
epochs l|Suntzefi. 1996. ). 

Starting from the photometric similarities of SN 2002dj 
and SN 2002bo, we studied the characteristics of other HVG 
SNe. Although the number of objects for which reliable 
multi-band photometry covering both very early and late 
light curves phases is limited, the HVG SNe analyzed in 
this paper seem to share the following photometric charac- 
teristics if compared with LVG SNe: 

1) A general fast rise to maximum brightness in all filters 

2) More pronounced inflections in the V and 7? bands 
around +25 days 

3) A brighter, more slowly declining B light curve after 
+40 days 

4) Possibly different colours 

To complete the picture these properties should be com- 
bined with parameters one can derive from the spectral evo- 
lution: 

1) Fast expanding ejecta 



2) Strong absorption due to intermediate mass elements 
(IME) at early phases 

3) Absence of carbon in the very early spectra 

4) Early time low photospheric temperature rapidly in- 
creasing toward maximum brightness 

5) Nebular spectra with kinematics not so different or 
even slower than that of LVG SNe. 

iBenetti et al.l l|2004h suggested that the high expansion 
velocities of IME are due to the fact that in HVG SNe the 
burning front extends further out into the outermost lay- 
ers than in LVG SNe. This hypothesis is supported by the 
extension of the blue wing of the Mg II line at velocities 
close to 19000 km sec~^ observable in the early IR spectra 
of SN 2002dj and by the absence of carbon in the -10.9 spec- 
trum of SN 2002dj and -12.9 spectrum of SN 2002bo. In this 
scenario the density of IME is enhanced in the layers where 
usually unburned carbon and oxygen are present. The IME 
density enhancement in the outermost layers would induce 
stronger absorptions at early epochs and increase the opac- 
ity in these regions, moving the photosphere to larger radii. 
The blending between high velocity features and the photo- 
spheric comp onent could also con tribute to the absorption 
enhancement (|Tanaka et al.|[2008l ). Since the photosphere is 
located at high velocity shells, this could justify both the fast 
rise to maximum observed in the light curve and, because 
the heating produced by the 7— rays is less effective, the 
lower temperature in the line forming region. Due to their 
high velocity the density in those shells decreases rapidly 
and hence the photosphere quickly recedes into the ejecta. 
Indeed, around maximum, the light curves of SN 2002dj, 
SN 2002bo and SN 2006X approached those of LVG SNe. 
Also the temperature is nearly identical as traced by the 
i?(Si II) evolution and by the appearance of the Si HI lines, 
the absorption minimum of which is at the same position in 
all SNe. The initially low and then rapidly increasing tem- 
perature of SN 2002dj was confirmed by fitting a black body 
to its two earliest IR spectra. 

We note that in HVG SNe more pronounced inflec- 
tions in the V and R light curves around +25 days seem 
to correlate with expansion velocities. The fastest expand- 
ing SN 1997bp, SN 2002bf and SN 2006X display a stronger 
inflection than in SN 2002bo and SN 2002dj, where the ef- 
fect is barely visible. We must stress that in the case of 
SN 2006X, due to the large reddening, the intensity of this 
feature has to be taken with caution. Also the brighter B 
light curve a fter +40 days seems t o correlate with the expan- 
sion velocity. IWang et all l|2008bl ) suggested that SN 2006X 
declined more slowly after +40 days and could be due to a 
light echo produced by the circumstellar dust. In the case 
of SN 2002dj there is no clear evidence of dust obscuring 
its light, but it is also the HVG SNe with the B light curve 
closer to "normal" SNe la. Nevertheless, SN 1997bp and es- 
pecially SN 2002bf with low absorption closely resemble the 
behavior of SN 2006X having also comparable expansion ve- 
locities. Moreover, part of the reddening estimated for those 
SNe could be biased due to their peculiar colour and colours 
evolution. 

Regarding the colours of HVG SNe we note that 
iDominguez. Hoflich fc Stranierol |2003) using a set of De- 
layed Detonation models (DD) have found that increasing 
the metallicity of the progenitor white dwarf, the B — V 
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colour of the resulting SN la becomes redder. Unfortu- 
nately, this study does not extend the results to other 
bands and is limited to a metallicity Z=0.02. However, 
iTimmes. Brown fc TruraiTI (|2003l ) have proposed that the 
effect of the progenitor metallicity is significantly enhanced 
at solar metallicity and above. In particular Fig. 9 of 
Dominguez et al. (2001) shows that passing from Z=0.001 
to Z=0.02 the flux in B is depressed, while that in V is en- 
hanced. If this is maintained at higher metallicity it could 
be the source of the red B — V and blue V — I colours we 
found in SN 2002bo, SN 1997bp and especially SN 2002dj 
where the effect of reddening should be negligible. Further- 
more, Timmes et al. (2003) have pointed out that in metal 
rich progenitors due to the more efficient electron capture, 
more ^*Fe and ^*Ni are produ ced at the expense of ^^Ni. 
iMazzali fc Podsiadlowskil (|2006l ) used this finding to explore 
the effect on spectra and light curves by varying the fraction 
of ^^Fe and ^*Ni produced in burning to nuclear statistical 
equilibrium and concluded that the colours of SN la with 
a larger fraction of these elements (i.e. the more metal rich 
ones) have redder B — V colour. Moreover, indications of 
larger '^'*Fe and ^**Ni and lower ''^Ni production for the HVG 
SN 1981B, SN 2002bo and SN 2002dj with respect the sim- 
ilar Ami5(B) LV G SN 1990N, SN 200 3 du and SN 2003cg 
have been found bv lMazzali et all (|2007l '). lLentz et all (|200ol i 
pointed out that increasing the abundance of all elements 
heavier than oxygen in the unburned C-l-0 layer of the 
pure deflagration model W7 (|Nomoto. Thielemann fc Yokoil 
1 19841 ') leads to a larger blueshift of SN la absorption lines. 
This is due to the fact that the higher opacity provided by 
the large amount of metals confines the line forming region 
to the outer (faster) part of the ejecta. Although, this sup- 
ports the hypothesis of metal rich progenitors for HVG SNe, 
we favour the configuration with an extended burning front 
to explain the HVG SNe observables. 

To further complicate the pi c ture w e must mention that 
iHoflich. Wheeler fc ThielemannI || 19981 ) found that SN la 
with metal rich progenitors should have a slower rise time 
to maximum than metal poor ones clearly contrary to what 
we found of SN 2002dj, SN 2002bo and SN 2006X. Given 
the very small sample of HVG SNe, the uncertainties af- 
flicting their colours and the model shortcomings we briefly 
recapped here, our suggestion of metal rich progenitor for 
HVG SNe are tentative and provide possible avenues to ex- 
plore when more data on HVG SNe and more sophisticated 
models will be available. 

F ast expanding SNe are rare both at low and high red- 
shift (|Blondin et all l2006l : iBalland et al] l2007l : iHook et all 
I2OO5I ). The systematic error they could introduce into cos- 
mological studies due to their peculiar colours appears to 
be negligible. Detailed analysis will become possible when a 
larger sample of well-observed rapidly evolving SNe la will 
be available. 
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Figure Al. Comparison of the different instrumental UBVRI 
transmission curves normalized to the peak transmission with the 
standard Johnson Cousins functions (,BessclL1990.V 



APPENDIX A: S-CORRECTION 

To avoid the well known limitations of the colour equation 
to calibrate SN, we used the following calibration path to 
transform the SN instrumental magnitude into a standard 
photometric system. To simplify the explanation of the var- 
ious steps, we will discuss as an example the calibration of 
the B band and we will use only one of the local standard 
stars. First one needs to transform the standard magnitude 
of the sequence star in the natural photometric system of 
the instrument using the following formula; 



(Al) 



where B* and V* are the standard magnitude of the local 
standard star, while 7b_v is the colour term. Then one com- 
putes the instrument zero point in the natural photometric 
system ZP{B)nat=B^at — b", where b* is the local standard 
star instrumental magnitude. 

Finally we calibrated the magnitude of the SN using the 
equation: 



= ZP{BUt + 6'''^ + 5b 



(A2) 



where b'^'^ is the instrumental magnitude of the SN. We did 
not correct for atmospheric absorption because we assume 
that the SN and local standard stars are at the same airmass 
and hence the absorption is already included in ZP{B)„at. 
We define as the difference between the synthetic magni- 
tude of the SN computed using the standard and the natural 
system passbands. 



Sb = Bt 



with: 



B 



-2.5log 



/_+^iV(A)gij(A)dA 



+ ZP{B)sy 



(A3) 



(A4) 



= -2.5^05 



/_+^jV(A)g^„,(A)dA 
/+^5f„,(A)dA 



-f ZP{h)s 



(A5) 



where A'^(A) is the photon number distribution of the 
SN, S^atW a-nd S'^(A) are the instrumental and standard B 



passbands. In the eauations lA4l and lA5l A'^(A) is used instead 
of -F'(A) because CC D and IR detector are photon counting 
devices. The BcssclJ l|l990t ) passbands were also reported in 
photon u nit dividing the orig inal ones for the wavelength, 
while the lPersson et ahl (|l998l ) bands are already in photon 
units. Since in the spectra which are used to compute the 
synthetic magnitudes the flux is expressed in energy units, 
7V(A) becomes: 



iV(A) = 



F{X)X 
he 



(A6) 



The zero points ZP{B)sy and ZP{b)sy were computed 
using a subset of spectrophotometric standard stars for 
which accurate photometry is available. In particular, for 
ZP{B)sy the standard magnitude was used directly, while 
for ZP{b)sy the magnitudes were previously reported in the 
instrument natural photometric system using equation HAlf) . 
Putting together equations lAll IA2I and IA3I we obtain: 



6*-7l-v(5'-^*)+S; 



'C(A7) 



Al Band construction 

We define the transmission function of a give passband 5(A) 
as: S(A) = F{X) ■ QE{X) ■ Ac{X) ■ Al{X) ■ M(A) • N^ef ■ L{X) 
where -F(A) is the filter transmission function, QE{X) is the 
detector quantum efficiency, Ac{X) is the continuum atmo- 
spheric transmission profile, Al{X) is the line atmospheric 
transmission profile, M(A) is the mirror reflectivity func- 
tion, Nref the number of reflections and L(A) is the lens 
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Table Al. Comparison between synthetic and photometric colour terms. In general we used B and B — V, V and B — V, R and 

V -Ras well as / and i? - /. But for EMMI and WFC we use V and V - R, for CTIO 0.9m and FORSl we use / and V - and 
for ALFOSC we use Rand R- I and I and V - I. 

sy^ ph^ sy ph sy ph sy ph sy ph 

CTIO 0.9m 0.062(0.035)" 0.098(0.013) -0.092(0.007) -0.094(0.007) 0.018(0.001) 0.016(0.003) 0.010(0.001) -0.013(0.009) 0.004(0.001) 0.004(0.004) 

ALFOSC 0.093(0.006)" 0.076(0.045) 0.032(0.003) 0.057(0.016) -0.046(0.002) -0.052(0.007) -0.065(0.002) -0.072(0.020) -0.041(0.002) -0.065(0.029) 

DFOSC 0.058(0.006)" 0.029(0.016) 0.068(0.003) 0.089(0.004) 0.013(0.002) 0.013(0.007) 0.026(0.001) 0.009(0.015) -0.048(0.001) -0.060(0.036) 

■IKT 0.072(0.050) 0.026(0.019) 0.066(0.006) 0.085(0.032) 0.030(0.001) 0.028(0.035) 0.003(0.001) -0.039(0.023) 0.056(0.004) 0.018(0.030) 

WFC 0.077(0.021) 0.059(0.032) 0.031(0.006) 0.063 -0.010(0.001) 0.007(0.029) 0.013(0.001) 0.000(0.014) -0.197(0.004) -0.210(0.004) 

FORSl - - -0.083(0.007) -0.083(0.009) 0.040(0.002) 0.033(0.010) 0.080(0.002) 0.066(0.010) -0.034(0.004) -0.031(0.005) 

EMMI - - 0.067(0.006) -0.045(0.005) 0.017(0.001) 0.056(0.023) 0.067(0.001) 0.044(0.009) -0.047(0.001) -0.038 

" Modified band. 

' Colour term computed using synthetic photometry. 
Colour term computed using photometric standard. 



Table A2. Comparison between synthetic and photometric colour terms for the infrared instruments. For SofI we use J and J — K, 
H and J - K, Ks and J - K, while for ANDICAM we use J and J - Ks, H and J - Ks, Ks and J - Ks. 



ay ph ay ph sy ph 



SoFi -0.020 -0.007 -0.022 0.006 0.006 0.023 

ANDICAM 0.035 0.028(0.005) -0.009 -0.010(0.005) 0.000 -0.003(0.005) 
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Figure A2. Comparison of the different instrumental JHKs 
transmission curves normali zed to the peak trans mission with the 
standard Persson functions llPersson et al.|[l998l) . 



a problem for the optical R and I filters, where some telluric 
features fall in the middle of the passbands. These telluric 
lines do not change linearly with airmass, therefore their 
variation is not properly accounted by the abs orption coef - 
ficients. For La Palma we obtained Ac{\) from lKind (llQSST ). 
while for La Silla, Paranal and CTIO we have used the CTIO 
transmission curve in the IRAF reduction package (Stone & 
Baldwin 1983, Baldwin & Stone 1984). For the IR sky spec- 
trum of La Silla and CTIO we downloaded Ac{X) ■ Al{X) 
from the 2MASS web sideQ. The instrumental filter bands 
are compared with the standard ones in Fig. lAll and Fig. lA2l 
As usual (see e.g. Pignata et al. 2004), we checked how well 
the constructed instrumental passbands match the real ones 
and compared the colour terms derived from standard stars 
(Landolt (1992) for the optical and Persson et al. (1998) for 
the IR) with those we compute d by integrating the fiu x of 
spectrophotometric standards of lStritzinger et al] (|2005ll for 
the optical and Sirius, Vega and the Sun for the IR. We de- 
cided to modify the instrumental passband only in cases 
where the difference between the two colour terms reported 
m Table and Table \M\ were larger than three times 
the error associated to both "synthetic" and "photometric" 
colour terms. The lack of a suitable number of spectrophoto- 
metric IR standards with an establish magnitude make the 
colour term comparison less effective than at optical wave- 
lengths. 



throughput. We downloaded F{X) and QE{X) from the in- 
strument web sites. For Ai (A) we used a standard aluminium 
reflectivity curve, while we did not find information about 
the lens transmission for any of the instruments. In the con- 
struction of the atmosphere model the weight of the terms 
Ac{X) and Al{X) change considerably going from U to Ks 
bands. In the optical the atmosphere transmission is mainly 
defined by a continuum function, while in the IR in most of 
the instruments the passband cut-offs are defined by deep 
atmospheric absorption bands. The latter make the use of 
a fixed instrumental passband not a good approximation of 
the real transmission function, because it depends on the 

atmospheric conditions. To a much lesser degree this is also * [ http://spider.ipac.caltech.edu/staff/waw/2iiiass/opt_cal/soutli_atm.tbl.htmII 
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Table A3. Optical photometry of SN 2002dj calibrated using the colour equation 
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" Relative to the time of the B maximum brightness M.J.D. =52450 ± 0.7 



Table A4. Near infrared photometry of SN 2002dj calibrated using the colour equa- 
tion. 



date 
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" Relative to the time of the B maximum brightness M.J.D. =52450 ± 0.7 



